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Abstract. In this paper, the dynamics of nonholonomic systems on Lie 
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are considered. Equations of motion form a closed system of differential 
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1 Euler — Poincare — Suslov nonholonomic systems 

Let the configuration space of a mechanical system with n degrees of freedom 
represent a Lie group G, and let wi,..., Wn be independent left-invariant vector 
fields on G (they are invariant under all left translations on G). The commutators 
of these fields have the form 


n 

[Wi,Wj] = 

/c=l 


( 1 ) 


where cf,- = —cj,- are the structure constants of G. 

Let X = {xi,... ,Xn) be local (generalized) coordinates on G and f{x) be 
a smooth function. Then 


/ 



n 

i=l 


where Wi (/) = -Wi is the derivative / along the vector field Wi. The variables 
oj (called quasi-velocities) linearly depend on x: 

n 

Xk = '^Wi{xk)uJi, l<k<n. (2) 

2 = 1 


The quasi-velocities are Cartesian coordinates on the Lie algebra g of G. 

Assume that the Lagrangian of the system of interest is reduced to the kinetic 
energy T which generates a left-invariant Riemannian metric (•, •) on G. In this 
case ^ ^ ^ 

T = - (x, i) = - ( ^ WiUi, '^ 3 ^ 3 ') = 2 X! 

where 

lij = {wi,Wj) = const, 


since all the vector fields ,..., Wn are left-invariant. The scalar-product matrix 
ll-^d II is positive dehnite. This is the inertia tensor of the mechanical system. 

According to Poincare [T^] (see also ID), the Lagrange equations take the 
following form: 

71 

Wi = ^ 1 <i <n. (3) 

3,k=l 


Here 


nii = 


dT 

duji 




are the components of the angular momentum of the system; this is a vector 
from the linear space g* dual to g. Eqs. ([3]) are a closed system of differential 
equations with quadratic right-hand sides on g (or on g*). The system is often 
called Euler - Poincare equations on a Lie algebra. For the group of rotations of 
the three-dimensional Euclidean space 50(3), Eqs. ([3|) coincide with the Euler 
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equations from rigid body dynamics. To completely describe the motion, one 
needs to add the kinematic equations Q to Eqs. ©• 

Now we increase the complexity of the problem by adding some constraints 
that are linear in the velocities and generally nonintegrable: 

fi{u!,x) = ... = fmioJ,x) = 0 (m < n). 


Of particular interest is the case where, apart from the kinetic energy T, the 
functions of the constraints are also left-invariant (i.e., explicitly 

independent of x). Under this condition the nonholonomic Lagrangian equations 
with the constraints 


ft fit 

c^rukUJj + X! 

j,k=l p=l ® 


= ... = fm{0j) = 0 


(4) 


form a closed system of differential equations on the Lie algebra g. 

For the case G = 50(3) systems with left-invariant constraints were first 
studied by Suslov [^. He considered the rotation of a rigid body about a fixed 
point with the following nonholonomic constraint: the projection of the angular 
velocity vector w € SO{3) = onto a body-fixed axis equals zero. General 
nonholonomic systems on Lie groups with a left-invariant kinetic energy and left- 
invariant constraints were studied in [201 E]. These systems are called Euler- 
Poincare - Suslov (EPS) systems. Since the constraints are linear in the velocity, 
Eqs. (HI) admit the energy integral 



E 




= h = const. 


(5) 


In Suslov’s problem Eqs. (H]) have the form 


lu! + u! X lu! = Aa, (a, ui) = 0. 


( 6 ) 


Here a ^ 0 is some constant vector in a moving space. If a is an eigenvector of 
the inertia operator {la = pa, p G M), then it follows from ([H]) that uj = const 
(see HS]). In particular, all rotations of the rigid body are stable. 

In a typical case (where this condition does not hold), each solution ti-^uj{t) 
possesses the following property: 

lim uj{t) = iwQ. (7) 


The motion of a rigid body is an asymptotic transition from steady rotation 
about some body-fixed axis to steady rotation with the same angular velocity 
about the same axis but in the opposite direction. We emphasize that in fixed 
space, as t —>■ -l-oo and t —oo, these rotations occur about different axes. The 
angular velocities of the steady rotations © are found from the system ([5]). On 
the one hand, we have the constraint equation (a, wq) = 0. On the other hand, 
performing scalar multiplication of the first equation ([6]) by /w°, we obtain 
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Fig. 1: 


one more equation: {Ia,uJo) = 0. This yields all steady-state solutions of the 
system m- 

u!q = fia X la, /r € K. (8) 

The phase portrait of the system is shown in Fig. [TJ 

It is evident from this phase portrait that the phase flow of this system 
does not admit an invariant measure with a continuous (and even summable) 
density. The conditions for the existence of an invariant measure of the EPS 
equations (|3]) in the multidimensional case are given in [SD]. Equality ([7]) holds 
also for Suslov’s problem in the n-dimensional space where G = SO{n) [TT] (see 
also nmn]). 

The absence of an invariant measure prevents one from applying the 
geometrical version of the Euler - Jacobi theorem which was formulated in [HIS], 
although explicit quadratures can be obtained (see [IS]). This motivated the 
development of an approach m based on the Hamiltonization of this system in 
open regions of the phase space. In [9| , it is shown that under certain conditions 
the complete system in Suslov’s problem (including the evolution of the Euler 
angles) can admit an additional algebraic integral of arbitrarily high degree in 
velocities. The explicit form of these integrals is found in where it is also 
noted that in this case the angle (in a fixed coordinate system) between the axes 
of limiting steady rotations is ±7r. In the general case, when the explicit form of 
first integrals is unknown, this issue is discussed in detail in m, and the angle 
between the same axes is determined by some formula depending on the system 
parameters. 

Using explicit formulas to solve Eqs. © |2S|, one can conclude that steady 
rotations of (I5|) are unstable for ^ > 0, and conversely, they are stable for ^ < 0 
(they are even asymptotically stable if the value of the total energy h > 0 
from ([S|) is fixed). 

Thus, the nonholonomic Suslov top can be spun about the axis (jSj) only in 
one direction: when spun in the opposite direction, the top loses stability and, 
with the course of time, begins to rotate in the opposite direction. This property 
of the Suslov top, usually called reversal, makes it similar to rattlebacks, which 
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also exhibit the asymmetry of stability when the spin direction is reversed m, 
see also [23]. We note that in [^ the property of reversal was noticed for the 
Chaplygin top (a dynamically asymmetric ball with a displaced center of mass 
rolling on a horizontal plane in a gravitational field). 

For recent developments in the study of rattleback dynamics, see IT6ll4l[^. 
where, in particular, it is shown that the absence of an invariant measure can 
lead to the existence of a system of strange attractors in phase space. 

The purpose of this note is to show the universality of the phenomenon 
of stability changes when the direction of the velocity of steady motions of 
nonholonomic EPS systems is reversed. Note that in the general case the absence 
or presence of tensor invariants mini El in nonholonomic systems leads to 
a kind of hierarchy of dynamical behavior described in 0 [30] • 


2 The Chaplygin sleigh as an EPS system 


The Chaplygin sleigh is a rigid body with a nonholonomic constraint moving on 
a horizontal plane: the velocity of some point of it is always orthogonal to the 
body-fixed direction (see [261 El])- For instance, one can assume that a vertical 
wheel which is unable to move in the direction orthogonal to its plane is rigidly 
attached to the body. We show that this nonholonomic system is an EPS system 
on the group of motions E(2) of the Euclidean plane. 

Without loss of generality, the rigid body itself may be assumed to be flat. 
The positions of this body on the plane are determined by three parameters: 
the Cartesian coordinates x and y of a distinguished point O, and the rotation 
angle (p (Fig. [2]). The independent left-invariant helds on E{2) are specified by 
the differential operators 


^ d . d 

= coup) — + sin(/j —, 
ox oy 


Xri = — sin ip 


dx 


cos if 


dy 
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The commutation relations ([T]) are 

[X^,X^]=-Xr,, [X^,X^]=X^, [X5,X^]=0. (9) 

Let Tj be the coordinates of the center of mass of the body G in a moving 
reference frame, u, v be the projections of the velocity of the origin of the 
moving frame onto these axes, and w be the angular velocity. The nonintegrable 
constraint equation is 

f = v = 0 . (10) 


This constraint is left-invariant. 

Let m be the mass of the body and Jq be its inertia moment relative to the 
center of mass. The kinetic energy 

is also left-invariant: like the constraint function, it is independent of the generalized 
coordinates. Hence, the Chaplygin sleigh is an example of an EPS system on 
the group E{2). 

Taking into account (jO]), the EPS equations for the Chaplygin sleigh take 
the form 


/dry dT /dry ar , 


ary ar ar 


( 12 ) 


Of course, one needs to add to them the constraint equation (nni). 

The second equation serves as the basis for finding the Lagrange multiplier 
(the constraint reaction), and the first and the third equations (taking into 
account the constraint) have the following explicit form: 

ii — rjuj = (Jo + rny -I- 1111 ^) 0 ; — mr]u = —m^uuj. (13) 


These equations should be supplemented with those describing the law of motion 
of the point O and the rotation angle ip: 

x = ucosip, y = usmip, ip = u!. (14) 

A qualitative analysis of motion of the Chaplygin sleigh for 77 = 0 can be found 
in |24| . We supplement this analysis with some remarks. 

Proposition 1. )1?]/ . Ifrj^O, then the equations of motion in the coordinate 
frame Oifrj with the origin 0\ lying at the intersection of the straight line OiG 
which is parallel to the blade and passes through the center of mass, and the 
straight line OOi which is perpendicular to OiG and passes through the contact 
point of the blade, coincide with the equations of motion (HI, (HI in the frame 
Ofrj when r] = 0. 
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Thus, if we consider the trajectory of point Oi instead of the trajectory of 
the contact point of the blade, then without loss of generality we may set 77 = 0 
in the equations of motion (ITSl) . 

If = 0, all motions are steady. If ^ 7 ^ 0, the phase portrait of the system (I13|) 
is similar in appearance to the phase portrait of Suslov’s problem (Fig. [T]). In 
particular, the nonlinear equations (1131) do not admit an invariant measure with 
a summable density. The steady-state solutions of m 

U = Uq, w = 0 


correspond to the rectilinear motions of the Chaplygin sleigh. They are stable 
for > 0 s-iid unstable for < 0. Thus, during a stable motion in a straight 
line, the center of mass of the body “outstrips” the contact point of the wheel. 
Here a change in stability occurs when the direction of motion is reversed. 

In fact, the system’s global evolution on the time interval t £ (— 00 , +c») in 
this case may be considered as a scattering process, that is, as t ^ — 00 , the 
sleigh “starts” its motion from some unstable steady-state solution (the center of 
mass “lags behind” the contact point), undergoes some evolution, and tends to 
some stable steady-state solution as f —?► -l-oo. At the same time, the complete 
rotation angle of the axis (or the axis O^Ci) is independent of the initial 
conditions and, according to HZ], is defined by 


= ttA, 

~ I —00 ’ 



The Hamiltonian representation of the system (HSl-dll and the redundant 
set of first integrals are also given in HZ). 

Typical trajectories of point Oi for various values of A are shown in Fig. [3l 
Examples of other nonholonomic systems represented in Hamiltonian and 
conformally Hamiltonian forms are given in [zi [3 la [3|. We also note that the 
motion of the Chaplygin sleigh on a rotating plane is considered in Ell- 


Remark 1. There is a simple and natural connection between Suslov’s and 
Chaplygin’s problems. We consider the generalized Chaplygin sleigh on a two- 
dimensional sphere: this is a spherical “cap” sliding on the sphere with the 
nonholonomic constraint. Its dynamics are described by Suslov’s equations. In 
the extreme case, as the radius of the sphere increases to infinity, we obtain 
the Chaplygin sleigh. This technique includes the well-known retraction of the 
group SO{S) to the group E{3). 


3 Anisotropic friction 

A nonholonomic model describes the dynamics of systems with ideal constraints. 
There exist more realistic models which take into account large forces of viscous 
friction with anisotropic Rayleigh’s function, see, for example, (TB]. The 
corresponding passage to the limit with regard to the Chaplygin sleigh was 
studied in [HI [21]. 
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Fig. 3: 


Following Caratheodory, we consider the motion of a rigid body on 
a horizontal plane taking into account the force of viscous friction applied to 
a fixed point of the rigid body and orthogonal to some body-fixed direction. In 
the notation of Section [2] the equations of motion are described by the following 
closed system on the algebra e(3): 


/ary _ ar 

\ au y ^ dv' 



dT , 

-UJ- - kv, 

ou 


dT\ dT dT 


(15) 


Here, k > 0 is the coefficient of viscous friction. In contrast to m, there are 
no constraints here. We write the explicit form of these nonlinear equations: 


u — rjuj = Luv + 

J 


V + = —oju + rjoj^ -u, 

m 


(16) 


— ui + — rjii = —uj{^u -I- Tyu). 

m 


Here, J = Jq + m(f^ -|- rf') is the moment of inertia of the body about the 
distinguished point. 

Like (jl2l) . the system m has the family of steady-state solutions 


u = Uq, V = 0, w = 0 


(17) 


corresponding to the motion of a rigid body in a straight line with a constant 
velocity. In view of the equality T = —kv'^, the energy does not dissipate on 
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these solutions. In order to study the stability of the solutions (113, we linearize 
in their neighborhood Eqs. dm). Omitting simple calculations, we give the 
characteristic polynomial of the linearized system: 


/(A) = A 


Jo . ^Jo + 


X^— + X 

m 


^uok 


(18) 


The appearance of the zero root is due to the fact that the steady-state solutions 
m are nonisolated. 

If ^uq < 0 (the contact point “outstrips” the center of mass when the body 
is moving), the characteristic polynomial has a real positive root and therefore 
the steady motion in a straight line is unstable (as in the case of the Chaplygin 
sleigh). 

We show that given > 0 (when the center of mass “outstrips” the contact 
point) the steady motion is stable. Moreover, perturbed motions asymptotically 
tend to one of the steady states (II3- By Lyapunov’s classical theorem 
Section 32], it is sufficient to show that the nonzero roots of the characteristic 
polynomial (1151) have negative real parts. But this, in turn, follows from the 
positivity property of the coefficients of the polynomial (1181) . 

We also present the asymptotic form of the roots of the characteristic 
polynomial as k ^oo: 


Ai = 0, A 2 


m^uo 
Jo + 



A 3 


Jo + 

Jom 


k + 0{l). 


It is clear that A 3 —> — 00 as fc —^ 00 , and 


A2 —t — 


m^uo 
Jo + ' 


(19) 


The unbounded decrease in one of the eigenvalues corresponds to the limit 
passage to the nonholonomic dynamics of the Chaplygin sleigh (taking into 
account the boundary layer phenomenon [ilM]). The limit relation (1191) gives 
a formula for the eigenvalue of steady-state solutions of the nonholonomic 
equations (fT3l) . 


4 The general case 

These observations can be generalized. Eqs. o are a closed system oi n — m = 
p > 2 differential equations with quadratic right-hand sides. Indeed, quasi¬ 
velocities can be chosen in such a way that the linear functions of constraints 
take the form: /i = ujn-m+i-, ■ ■ ■, fm = Then the first n — m equations of (U) 
form a closed system with the quadratic right-hand side 

UJ = v{uj), UJ = ( 20 ) 

and Lagrange multipliers can be found from the remaining m equations. The 
system (1501) admits the first integral iL as a positive definite quadratic form. 
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The problem of finding nonzero steady-state solutions ujq to the nonlinear 
system (l20l) is a nontrivial algebraic problem. It depends on the structure 
constants of the algebra g, the inertia tensor, and the linear functions of 
constraints. The homogeneity of the system (uni) 

v(Auj) = A^v(uj), A € M, (21) 

implies the following simple property: if a; = wq is a steady-state solution, then 
the whole straight line 

w = g,ujQ, /i € K., (22) 

consists of steady-state solutions. We call the straight line (|2^ a straight line 
of steady motions (which we abbreviate as SLSM). 

The general number of nonzero nonproportional solutions of the algebraic 
system of equations 

ui(a;) = 0,..., Vp(w) = 0 (23) 

(including the complex solutions) can be estimated using the Bezout theorem: 
if the number of these solutions is finite, it does not exceed 2^“^. Indeed, from 
the algebraic system (1^51) we select p — 1 equations so that the number of their 
nonproportional solutions (with the complex ones) is finite. Since the number 
of variables equals p and all functions vi,... ,Vp are homogeneous in w with 
homogeneity degree 2, then (by the Bezout theorem) the number of solutions 
(with multiplicities taken into account) equals 2p~^. All these solutions should 
satisfy one more (omitted) equation. Therefore, their general number obviously 
does not increase. 

This estimate is rough and does not take into account the existence of the 
positive definite quadratic integral of the system (I20II . 

Theorem 1. Ifp = 2, then either all positions are steady {vi = U 2 = 0) or there 
is only one SLSM. 

Thus, if p = 2, the phase portrait of the system ([20ll is similar in appearance 
to the phase portraits in Suslov’s and Chaplygin’s problems. Theorem [T] shows 
universality of stability changes in EPS systems with two nonholonomic degrees 
of freedom when the direction of motion is reversed. Theorem [T] also shows that 
the general estimate of the number of SLSMs based on the Bezout theorem is 
not accurate. 

Proof, of Theorem [TJ By a linear change of independent variables the first 
integral of the system (EOl) is reduced to the form 

^ = 2 (*^1 ^ 2 )- 


In these variables, let 


Vi = aiUj\ + biU}iU}2 + CiUi\, i = 1, 2. 
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The condition H = 0 gives the following relations for the coefficients: ai = C 2 = 
0, bi = — 02 , Cl = — 62 . Hence, Eqs. (l20l) become 


OJl = UJ2{blUJi + C1UJ2), W2 = —UJi{biUJi + CiW2). 


(24) 


ff bi = Cl = 0, then vi = V 2 = 0. But if bl + cl ^ 0, there is only one SLSM: 
biUJi + C 1 W 2 = 0. QED. 

Eqs. (1^^ can be given the conformally Hamiltonian form 


oil 



102 = -P 


dH 

duji ’ 


p = biUJi + CiUl 2 - 


They take the canonical form after rescaling time asdr = pdt. The multiplier p~^ 
serves as the density of the integral invariant of the system (1^01) . However, this 
function has singularities on the SLSM and reverses sign when crossing this 
straight line. Conformally Hamiltonian systems naturally arise in problems of 
nonholonomic mechanics [21 [3]- 

Let w = wq be a steady motion. Set 


A = 



(25) 


This is the linearization operator of the system (EUl) at point loq. Differentiating 
the identity (1211) with respect to A and then assuming that A = 1, we obtain 


If w = tuo, we have Aloq = 0. Therefore, the operator (1^51) is degenerate. By ^(wo) 
(s(coo)) we denote the number of eigenvalues of the operator (1251) lying in the 
right (respectively, left) complex half-plane. The following simple theorem holds: 

Theorem 2. If coo ^ 0 is a steady motion, then for /i > 0 


u{pLOo) = u{loo), s{pLOo) = s(wo), 

and for p < 0 

u{ploq) = s(wo), s{pujo) = u{ujo)- 

Indeed, according to (I25L the eigenvalues of the linearization operator of at 
the stationary point poog differ from those of the operator A by the multiplier p. 
In the case that is the most important to us 


s(u;o) = n — m — 1. 


Then (with the zero root) u(coo) = 0 and, by Lyapunov’s theorem, the equilibrium 
point loq is stable. Moreover, the perturbed motion infinitely approaches one of 
the equilibria (l22l) as t — )• -|-c». On the other hand, by Theorem |2l tt(— wq) also 
equals n — m — 1 > 0. Therefore, the steady motion co = —loq is unstable. 
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We call the SLSM nondegenerate if at all points of this straight line, 
oj ^ 0, the operator ((25l) has only one zero eigenvalue. We define the index of 
the nondegenerate SLSM: 

* = j(wo) + Wo 7^0, (26) 

where j{u!o) = 1 (= —1) if s(u;o) is even (odd). By Theorem [21 the index is 
independent of the choice of point uo ^ 0 on the straight line (l22l) . The index 
can take the following values: —2, 0, 2. If p is even, then (by Theorem|2]) i = 0. 

Theorem 3. If all SLSMs are nondegenerate, then the sum of their indices 
equals 1 + (—1)^+^. 

In order to prove it, we fix the positive level set of the energy integral 
H = h > 0 and restrict the initial dynamical system (1201) to it. It is only the 
equilibrium points of the reduced system on the {p — l)-dimensional ellipsoid 
{Hiuj) = h} that coincide with the points of intersection of the SLSM with 
this ellipsoid. In view of the assumption about nondegeneracy of the SLSM, all 
equilibria of the reduced system are isolated. The index of the singular point wq 
of the system on the energy ellipsoid is defined as the sign of product of the roots 
of the characteristic polynomial of the linearization operator for the reduced 
system at point wq. It is clear that these numbers coincide with the eigenvalues 
of the operator (1^51) other than zero. Their product sign obviously coincides 
with the number j{uJo) from (1261) . It remains to use the Poincare-Hopf theorem: 
the sum of the indices of isolated singular points of a dynamical system on 
a closed manifold is equal to its Eulerian characteristic. Recall that the Eulerian 
characteristic of the {p — l)-dimensional sphere is equal to zero for even p and 
is equal to 2 for odd p. QED. 

Corollary 1. If p is odd, then there is at least one SLSM. 

Corollary 2. If p is odd and all SLSMs are nondegenerate, then their number 
is odd. In particular, it does not exceed 2P~^ — 1. 

Indeed, for odd p the index of the SLSM can be either —2 or 2. If the number 
of nondegenerate SLSMs is even, then the sum of the indices is divisible by 4. 
However, (by Theorem [3]) this sum equals 2. Further, by the Bezout theorem, 
there are at most 2p~^ nondegenerate SLSMs. But for p > 2 this is an even 
number. Hence, the total number of SLSMs does not exceed 2^“^ — 1. 

It is not improbable that the conclusions of Corollaries |T] and |2| hold for even 
values of p, too. In any case, this holds for p = 2 (Theorem [1]) . When p = 3, 
there is either one or three nondegenerate SLSMs. The latter case takes place 
for the Euler top when the inertia moments are unequal. 

In conclusion, we mention that a nonholonomic system of two coupled bodies 
(called a nonholonomic hinge and generalizing the Suslov problem) is considered 
in [321 • Evidently, the above methods can be applied to this problem as well. 
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